Optimal quantum measurements for spin-1 and 

spin- 1 particles 



P.K.Aravind 

Physics Department, Worcester Polytechnic Institute, Worcester, MA 01609, USA 

February 1, 2008 

Abstract 

Positive operator valued measures (POVMs) are presented that allow 
an unknown pure state of a spin-1 particle to be determined with optimal 
fidelity when 2 to 5 copies of that state are available. Optimal POVMs 
are also presented for a spin-| particle when 2 or 3 copies of the unknown 
state are available. Although these POVMs are optimal they are not 
always minimal, indicating that there is room for improvement. 

1 Introduction 

The problem of determining an unknown quantum state of which only a finite 
number of copies is available is one that has attracted much attention in re- 
cent years. A fruitful method of attacking the problem was suggested by Peres 
and Wootters[l], who pointed out that a judicious joint measurement on all 
the copies can sometimes yield more information than any measurements on 
the individual copies. This suggestion, whose formal basis is embodied in the 
theory of Positive Operator Valued Measurements (POVMs) [2], has served as 
the springboard for almost all subsequent work on the problem. Initial work 
focussed on determining an unknown (pure) state of a qubit of which N copies 
were available [3- 5], but this work was later extended to mixed states of qubits[6] 
and to pure states of qudits (i.e. D-state quantum systems) as well[7,8]. The 
relationship of quantum state estimation to quantum cloning machines was dis- 
cussed in [9]. Other aspects of the state estimation problem have also been 
discussed recently [10]. 

The purpose of this paper is to present some new POVMs that allow opti- 
mal determinations of unknown qudit states in certain cases. Section II presents 
POVMs that allow an unknown pure state of a spin-1 particle to be determined 
with optimal fidelity when 2 to 5 copies of that state are available. Section III 
presents optimal POVMs for a spin-| particle when 2 or 3 copies of the unknown 
(pure) state are available. The POVMs for spin-1 are based on the geometry of 
the 24-cell and the 600-cell, two four-dimensional regular polytopes, while the 
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POVMs for spin-| are based on the 40 states of the "Penrose dodecahedron" [11] 
as well as another configuration of 60 states. The geometric structures under- 
lying the various POVMs presented here have been exploited earlier to provide 
proofs of the Bell-Kochen-Specker[12] and Bell theorcms[13,14], and it is inter- 
esting that they should prove of use in the state estimation problem as well. 

2 POVMs for spin-1 

Given N identical copies of a Z?-state quantum system in an unknown pure state, 
how can one determine that state as reliably as possible? We first address this 
problem for general N and D, introducing our terminology and notation in the 
process, before specializing to the case D = 3 in this section and D = 4 in the 
next. The (now) standard method of attacking this problem is to implement a 
POVM (or generalized measurement) on the N copies and to use the outcome of 
this measurement to make a judicious guess about the input state. The success 
of a guess is gauged by a " fidelity" , which is generally taken to be the squared 
modulus of the overlap between the input state and the guess made for it. The 
problem of state estimation consists of devising a POVM for which the average 
fidelity (i.e. the fidelity averaged over all possible occurences of the input state) 
is as large as possible. If the input state is distributed with uniform probability 
over the entire projective Hilbert space of the qudit, it has been shown[7,8] that 
the average fidelity is bounded from above by the quantity ^xjj ■ The task of 
constructing a POVM that achieves this upper bound for any N and D was 
addressed in fairly general terms in ref. [8] , with the analysis for the cases N = 2 
and 3 (but arbitrary D) being carried somewhat further. However only for the 
case D = 3 and TV = 2 was a POVM explicitly constructed. This section will 
present optimal POVMs for D = 3 and N = 2- 5. 

An optimal POVM[2] for arbitrary N and D is characterized by a set of 
positive numbers c r and system states |\& r ) (r = 1, .., k) such that 

J2c r \* r f N N *{* r \=I, (1) 
r=l 

where |\l' r ) <g,Ar denotes the tensor product of the state |^ r ) with itself N times 
and I is the identity operator in the maximally symmetric subspace of the space 
of N qudits. The number of states in the POVM (= k) can be larger than the 
dimensionality of the maximally symmetric subspace (= ^tt^tw ) m which 
the POVM acts. A standard von Neumann measurement can be regarded as a 
POVM in which all the c r are unity and the operators that effect the resolution 
of the identity are projectors onto non-overlapping subspaces. A POVM can 
be implemented in practice by coupling the system of interest to an auxiliary 
system (the "ancilla") and carrying out a von Neumann measurement on the 
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enlarged system; this measurement then appears in the space of the system 
alone as a POVM[2]. 

The operator identity (1) can be turned into a scalar equation by taking its 
expectation value in )® N , the iV-fold tensor product of the arbitrary system 
state I* ). One then finds that 

Y^cr N9 (*\* r ) 9N N ®(y r \y)® N = j2 Cr (k*i^)i 2 ) jv = i - ( 2 ) 

r—l r—1 

We now specialize the discussion to a spin-1 system for which the spin value, J, 
is equal to 1 and the Hilbert space dimension is D = 2 J +1 = 3. The projective 
Hilbert space of a spin-1 system has dimension 2D — 2 = 4 and an arbitrary 
state, ), in this space can be parametrized as[15] 

|*) = (e lXl sin6>cos</>, e % sin6»sin0, cos6») where 0<Xi,X2< 27r > < 9, (f> < ir /2. 

(3) 

If one introduces the real variables x\ — cos <f) cos X\ 7 x i = cos <t> sm Xi 7 ^3 = 
sin <f> cos Xi an d x 4 = sm <P sm X2 that define a point on the surface of the four- 
sphere x\ + x\ + x\ + x\ = 1, the state |\& ) can be expressed as 

\ty) = (sin# (xi + ix 2 ) , sin# (x 3 + ix 4 ) , cos6>) , < 9 < tt/2, — 1 < x\, X2, X3, X4 < 1. 

(4) 

Our strategy for constructing an optimal POVM for D = 3 is as follows. 
For each c r and 9 r we choose a set of points xu, x 2 i, x^i, xu [i = 1, .., n) on the 
surface of the four-dimensional unit sphere (the same set turns out to suffice 
for every r) in such a way that when the sum in (2) is carried out over all i 
for a fixed r, the dependence on the angular variables <f>, \\ an d Xi cancels out, 
leaving a function of 9 alone. The sum on the left side of (2) then reduces to 
a polynomial of degree k — 1 in cos 2 # that depends parametrically on the as 
yet undetermined elements c r and 9 r (r = 1, .., k) of the POVM. We finally nail 
down the POVM by choosing k arbitrary angles 9 r and fixing the constants c r 
in such a way that the polynomial on the left side of (2) reduces identically to 
the unity. 

It remains for us to specify how the points xu, x 2 i, x^, x^ (i = 1, .., n) should 
be chosen so that the cancellation over the angular variables <p, \\ an d X2 can 
be accomplished on the left side of (2). For V = 2 or 3 it turns out that 
choosing the 24 vertices of a 24-cell does the trick, while for V = 2 — 5 choosing 
the 120 vertices of a 600-cell does the trick (the 24-cell and 600-cell are two 
of the four-dimensional regular polytopcs[16]). The unit vectors to the vertices 
of a 24-cell are the 16 vectors |(±1, ±1, ±1, ±1) and the 8 vectors (±1,0,0,0), 
(0, ±1,0,0), (0,0, ±1,0) and (0, 0,0, ±1) (i.e. the vertices of a hypercube plus 
those of a hyperoctahedron), while the unit vectors to the vertices of a 600-cell 
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are the 24 vectors just mentioned plus the 96 vectors obtained by taking all even 
permutations of ^ (±r, ±1, ±r _1 , 0), where r = ^(1 + ^/5) is the golden mean. 

We now expand upon the above explanation to show how our optimal POVMs 
for N = 2 — 5 are constructed. The POVM for TV = 2 is based on three numbers 
(or "weights") ci, c 2 and C3 and three angles 6\, 9 2 and 6> 3 , with each angle giv- 
ing rise to 24 states based on the geometry of a 24-cell as indicated in eqn.(4). 
On using this POVM in the left side of (2) and taking 6>i = 7r/4, 9 2 — tt/2 and 
#3 = 0, we find that the summation can be carried out analytically to yield the 
polynomial expression 

2ci + 8c 2 + (8ci - 16c 2 ) cos 2 9 + (-4ci + 8c 2 + c 3 ) cos 4 6>. (5) 

For this to be identically equal to 1, it is necessary that the constant term be 
equal to 1 and that the coefficients of cos 2 8 and cos 4 9 both vanish. This is 
achieved by the choice of weights c\ — i, c 2 = yj and C3 = 0. We have therefore 
constructed an optimal POVM for N = 2 which is characterized by two weights 
ci = i and c 2 = yj and a total of 48 states (24 arising from each of the angles 
9 1 = tt/4 and 9 2 = tt/2). An optimal POVM for N ~ 3 can be constructed in a 
similar manner. This POVM involves four weights Ci,c 2 , 03,04 and four angles 
$1j$2j$3> $4 and contains 4 x 24 = 96 states obtained by combining each 9 
value with the vertices of a 24-cell in the manner indicated in (4) . On using this 
POVM with 6»i = tt/6, 6 2 = tt/4, 6> 3 = tt/3 and 64 = tt/2 in the left side of (2) 
we obtain a third-degree polynomial in cos 2 9 which reduces identically to 1 if 
one chooses the weights as Ci = J^, c 2 = yg, C3 = | and C4 = y|g . Our optimal 
POVM for N — 3 is therefore characterized by 4 weights and 96 states. 

The above construction cannot be used to obtain POVMs for N > 4, because 
the use of a 24-cell no longer allows the angular dependence on <p, \\ an d X 2 to 
be cancelled out on the left side of (2). However, as mentioned earlier, replacing 
the 24-cell by a 600-cell allows this cancellation to be accomplished for all N 
values from 2 to 5. For N = 2 and 3 the 600-cell yields POVMs that are not 
as economical as the ones found earlier, so we skip over these cases and pass 
to N = 4 and 5. The construction of the POVMs for N — 4, 5 proceeds in the 
same manner as explained earlier for N = 2 and 3, but with the difference that 
the role of the 24-cell is now taken over by the 600-cell. Each 9 value in the 
POVM gives rise to 120 states, and the total number of states in the POVM is a 
multiple of 120. The details of our POVMs for N = 4 and 5 (i.e. the weights c, 
and angles 9i) are summarized in the bottom half of Table I, while the top half 
summarizes our results for N — 2 and 3. For N — 2 ref.[8] presents a POVM 
with just 8 states, which is more economical than the one found here. However 
for N > 3 no optimal POVMs of any kind have been reported earlier. 

3 POVMs for spin-3/2 

Before presenting our POVMs for spin-|, we return to eqn.(2) and recast it 
in an alternative form based on the Bloch vector description of qudit systems. 
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In this description [15], the state vector |^ ) of a D-state quantum system is 
represented by a unit vector - the generalized "Bloch" vector n - in a real 
(D 2 — 1) -dimensional space. The inner product of two state vectors is related 
to the scalar product of the corresponding Bloch vectors by the equation 

K*il*i>l 2 = 2jVi (1 + 2Jfii " (6) 

where the spin value J is related to D by D — 2 J + 1. On substituting (6) 
into (2) and performing some manipulations, (2) can be recast in the form of a 
heirarchy of equations that must be satisfied by the POVM elements (c r , n r ) 
for r = 1, .., k. We now turn to these equations[8]. 

For arbitrary D = 2 J + 1 and N = 2, the equations to be satisfied by the 
POVM elements are[8] 



^c r = (2J+l)(J + l), (7) 

r=l 



^2 c, r (n r • n s ) = 0, (8) 



k ^ 2 J -|- l 

and ^c r (n r - n s ) 2 = ^ . (9) 

r— 1 

The two last equations must be satisfied for each value of s, so the above re- 
quirements amount to a total of (2k + 1) equations. 

For arbitrary D = 2 J + 1 and N = 3, the equations to be satisfied by the 
POVM elements are [8] 



^ Cr= (2J + 3)(2J+l)(J+l) ! (iQ) 



J^c r (n r • n s ) = 0, (11) 

r=l 



± cASr . R ,f=2£±mi±v, (12) 

r=l 
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a ST - ^ (2J+1)(2J-1) 
and 2^c r (n r - n s ) = ^ ■ ( 13 ) 

r=l 

Again the last three equations have to be satisfied for each s, leading to a total 
of (3fc + 1) equations altogether. It is interesting to note, from a comparison 
of (7)-(9) with (10)-(13), that any POVM for three copies can be turned into 
one for two copies simply by reducing each of the constants c r by the factor 
(2 J + 3)/3. This observation will prove to be of use below. 

We now introduce the 40 states of the "Penrose dodecahedron" [11] and show 
how they can be used to construct POVMs for a spin-| particle. Twenty of these 
states (termed "explicit rays" in [11]) are the spin +5 projections of a spin-| 
particle along the twenty directions from the center of a regular dodecahedron 
to its vertices. The remaining twenty states (termed "implicit" rays in [11]) are 
also associated with the vertices of the dodecahedron and have the property 
that the implicit ray associated with any vertex, together with the explicit rays 
associated with the three neighboring vertices, constitute a mutually orthogonal 
set of states. Explicit expressions for all 40 Penrose states are given in Table II, 
in the basis afforded by four of these states. From this table it is easily verified 
that each state is orthogonal to exactly 12 others and makes the same, constant 
angle with the 27 states it is not orthogonal to. This implies that the scalar 
products of Bloch vectors for pairs of Penrose rays have only the two values 

n r • n s = — for orthogonal rays, (14) 



or n r • n s = - for non-orthogonal rays. (15) 

Using (14) and (15) it is readily verified that a POVM satisfying (10)-(13) with 
J = I is obtained by taking all forty Penrose rays with a common weighting 
factor of c r = | for each. Using the remark just below (13), it follows that 
a POVM satisfying (7)-(9) is obtained by taking all forty Penrose rays with a 
common weighting factor of c r = \ for each. The correctness of these POVMs 
can also be checked directly, but much more tediously, by verifying that they 
make the left side of (2) reduce identically to the unity for an arbitrary choice 
of system state |\&) . 

We have discovered yet another POVM for J = | and N — 2 or 3, consisting 
of the 60 states in Table III. For N = 2 these states are to be taken with the 
common weighting factor of c r = |, while for N — 3 they are to be taken with 
the common weighting factor of c r = -|. The correctness of these POVMs can 
be verified in two distinct ways. The first is to note that each of these 60 states 
is orthogonal to 15 others, makes a constant angle with 32 others, and a second 
constant angle with the remainingl2 others. In other words, if ni denotes the 
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Bloch vector of any one of these states and n r ranges over the Bloch vectors 
of all the others, the scalar product of ni with n r can assume only one of the 
following three values: 



1 

ni • n r = 



for 15 orthogonal rays, 



(16) 



ni • ri r = 



for 32 non-orthogonal rays 



(17) 



and n! • n r = - 



for 12 non-orthogonal rays. 



(18) 



Using these scalar products in (7)-(9) and (10)-(13) allows us to confirm that 
the 60 rays, taken with the appropriate weighting factors, are indeed POVMs 
for N = 2,3. The second way of checking the correctness of these POVMs is 
to substitute them into the left side of (2) and show that it reduces identically 
to the unity for an arbitrary choice of system state |^). We say a word about 
the origin of this POVM: the first 24 states in Table III were introduced by 
Peres[17] and used by him to prove the Bell-Kochen-Specker theorem; we added 
the remaining 36 states to obtain a set satisfying all the conditions for a POVM. 

This completes our presentation of optimal POVMs for D = 4 (or J = |) 
and N — 2,3. We have not succeeded in finding any optimal POVMs for D = 4 
and N > 4. 

For N = 2 the Penrose rays provide an optimal, but perhaps not minimal, 
POVM since Acin et.al.[8] constructed a smaller POVM involving only 15 Bloch 
vectors equally inclined to each other; however Acin ct.al. did not demonstrate 
that physical states corresponding to these Bloch vectors exist, so their solution 
cannot be regarded as complete. For N = 3 the Penrose rays definitely provide 
a minimal POVM because it was shown by Acin et.al.[8] that a minimal POVM 
in this case cannot consist of fewer than 40 states. Our 60-state POVM, though 
not minimal, is nevertheless of interest because of the neat solution it provides 
to the same problem. 

4 Some open questions 

The optimal POVMs for spin-1 presented in this paper all involve a large num- 
ber of states. For N = 3, for example, our POVM involves 96 states whereas 
the minimal number is expected[8] to be a little above 18. We suspect, there- 
fore, that there is considerable room for improvement in our POVMs, as far as 
their economy is concerned. Realizing this improvement poses an interesting 
mathematical challenge, but also one that is not without physical interest since 
more economical POVMs would likely lead to more streamlined experiments. 
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The POVMs presented in this paper are "special purpose" ones, tailored 
to specific spin values and small numbers of copies. It is clearly desirable to 
generalize the algorithms presented here and devise POVMs for any spin value 
and any finite number of copies. In the case of qubits such a general algorithm 
was proposed in ref. [4] , where it was shown how a POVM for N copies can always 
be constructed out of a suitable set of (N+ 1) 2 states. The technique underlying 
this construction is to distribute unit vectors over the Bloch sphere in such a 
way that a cancellation over the azimuthal angle <f> is first achieved, following 
which a cancellation over the polar angle 9 is achieved by suitably adjusting 
the weights, c r , in the POVM. A very similar technique of angular cancellation 
(first over the "azimuthal" angles <f>,Xi an d Xi-> followed by a cancellation over 
the "polar" angle 9) was employed in the construction of our special purpose 
POVMs for spin-1. We suspect that this technique can be generalized to yield 
POVMs for any spin and any number of copies, but the specific way to do this 
has so far eluded us. 

The problem of realizing the optimal POVMs proposed here as von Neumann 
measurements on an enlarged space of the system and an ancilla is an interesting 
one worth addressing. This would bring the scheme of optimal measurements 
one step closer to experimental realization and also help highlight any problems 
connected with its practical implementation. 

Acknowledgement. I would like to thank Rolf Tarrach for stimulating 
my interest in the problem of optimal quantum measurements and for several 
valuable discussions on the subject. 
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N 


Angles r 


Weights c r 


# states 


2 


0i = tt/4, 2 = tt/2 


2 Cl ^' C2 = 2^ 7 1 


48 


3 


9 1 = tt/6, 2 = tt/4, 3 = tt/3, 4 = 


c l 27 ' C2 18 ' C 3 9 ' C4 108 


96 


4 


01 = tt/6, 2 = tt/4, 3 = tt/3, 4 = tt/2 


Cl = 45 ' C2 — ^' C 3 — T5 ' C4 = 3fi0 


480 


5 


0! = tt/6, 2 = tt/4, 3 = tt/3, 4 = tt/2, 
5 = tt/8,06 = 3tt/8 


Cl — 225 ' C 2 — 300 ' c 3 — 75 i c 4 — Jgoo ' 
_ 2-v^ _ 2+^2 
L5 - fin , o fin 


720 



TABLE I. Optimal POVMs for estimating an unknown (pure) 
state of a spin-1 particle of which N = 2 —5 copies are available. The 
POVMs for N = 2, 3 are based on the geometry of a 24-cell, while those 
for N = 4, 5 are based on the geometry of a 600-cell. The states of the 
POVMs are constructed by combining the values in column 2 with 
the vertices of a 24- or 600-cell in the manner indicated in eqn.(4), the 
total number of states in the POVM being either 24 or 120 times the 
number of values in column 2 (this total is indicated in column 4). 
The "weights" associated with all states sharing a common value 
are indicated in column 3. 
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<^} = (l,p,p 2 ,0) 


= (1,0,0,0) 


|* B ) - (0,1,0,0) 


|* B ) = (0,0,1,0) 




* L ) = (-1,0,^,1) 


|* G ) = (0,-l,p,l) 


|*c> = (P", 1,0,1) 


|*r>) = (p, 0,1,1) 




vj) = (0,P , 1,-1) 


= (l,p- 2 ,0,l) 


= (0,p,-l,l) 


|*m) = (p" 1 , 0,1,1) 




**) = (l,0,p,-l) 


|*/) = (l,p 2 ,0,l) 


|* P ) = (p- 2 , 1,0,1) 


|* Q ) = (0,l,p 2 ,-l) 




*s> = (1,1,0,1) 


|*jv) = (0,1,1,-1) 


= (-1,0,1,1) 


|*r> = (P 2 ,P,1,0) 




*' A ) = (0,0,0,1) 




* F \ = (0,p 2 ,l,p) 




= (p,0,l,p 2 ) 




* £ > = (p- 2 ,p 2 ,o,i) 




*^ = (0,l,l,p) 




* G ) = (i,o,-i, P - 1 ) 




*c> = (1-0, -1,P) 




* D ) = (l,l,0,p 2 ) 




*' 7 } = (l,l,0,p- 2 ) 




*^) = (0,l,l,p- 1 ) 




*«> = (-i,i,jr\o) 




* m ) = (-i,i,p,o) 




* ff ) = (p 2 , -1,1,0) 




*j) = (p, 1,-1,0) 








*n) = (l,P,l,0) 




*s) = (l,P 2 ,0,p- 2 ) 




*^) = (0,l,p 2 ,p) 




^)-(p,o,p 2 ,i) 




*t) = (1,-1,1,0) 



TABLE II. The 40 states of the "Penrose dodecahedron". The ket 
vectors corresponding to the states are shown as row (rather than 
column) vectors for convenience. The first 20 states, without primes, 
are the "explicit rays" while the last 20 states, with primes, are the 
"implicit rays". Each state bears a subscript indicating the vertex 
of the dodecahedron with which it is associated (see the papers in 
ref.ll for a picture of a dodecahedron with its vertices labelled with 
the letters used here). The states \^f) , \^e) , \^b) an d , consist- 

ing of an implicit ray and the explicit rays at the three neighboring 
dodecahedron vertices, form a mutually orthogonal set and are used 
as a basis in which the components of all the states are expressed. 
The symbol p stands for exp(i^) and a normalization factor of 1/V3 is 
omitted from many of the states. Note that each state is orthogonal 
to exactly 12 others and makes a constant angle with the remaining 
27. 
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*i) = (1,0,0,0) 




* 2 ) = (0,1,0,0) 


|* 3 ) = (0,0,1,0) 




* 4 ) = (0,0,0,1) 




* B > = (1,1,1,1) 




*6> = (-1,1,-1,1) 


I*?) = (-1,-1,1,1) 




*8> = (1,-1,-1,1) 




*9> = (1,1,1,-1) 




*io) = (1,-1,-1,-1) 


|*n) = (1,-1,1,1) 




tfia) = (1,1,-1,1) 




*i 3 ) = (1,0,1,0) 




*u) = (0,1,0,1) 


l^is) = (1,0,-1,0) 




*i 6 ) = (0,1,0,-1) 




*i 7 ) = (1,1,0,0) 




*i 8 ) = (1,-1,0,0) 


|*1 9 ) = (0,0,1,1) 




* 20 ) = (0,0,1,-1) 




*2i) = (-1,0,0,-1) 




*22> = (0,-1,-1,0) 


|*23> = (-1,0,0,1) 




*24> = (0,-1,1,0) 








*26> = (l,-»,-*,l) 


1*27) = (1, — — 1) 




*28> = (1,0-0-1) 




*2 9 ) = (— 1, 1, — i, — «) 




*3o> = (-1,-1, 


1*31) = (-1,-1, -00 




*32) = (-1,1,00 




*33> = (-1,-1,1,-*) 




*34) = (— 1,*, — 1,— *) 


1*35) = (— 1,— *,— 1,*) 




*36> = (-1,0 00 




*37> - (1,0,0,Z) 




*38> = (1,0,0,-i) 


|*39) - (0,1, i, 0) 




*4o) = (0,1,-00) 




*41> = (1,0, t,0) 




*42> = (1,0, -*,0) 


|*43) = (0,1,0,0 




* 44 ) = (0,1,0,-0 




*45) = (1,00-1) 




*4 6 > = (l,-i,-i,-l) 


1*47) = (1,0-01) 




*48) = (1,-001) 




* 49 ) = (l,i,0,0) 




*so) = (l,-i,0,0) 


|*5i) = (0,0,1,0 




* 52 ) = (0,0,1,-0 




* 53 ) = (1,*,1,») 




* 54 ) = (1, -i, 1, -i) 


1*55) = (1,0-0-0 




*56> = (1,-0-00 




* 5 ?) = (M,00 




*58> = (1,-1, 0~0 


1*59) = (1,1,-0-0 




*6o) = (1,-1,-0 



Table III. A set of 60 states yielding optimal POVMs for J = § 
and N = 2,3 (note that most of these states are unnormalized). The 
POVM for N — 2 is obtained by taking these states with an equal 
weighting factor of c r = while the POVM for N = 3 is obtained by 
taking these states with an equal weighting factor of c r = |. 
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